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On the rotation of ONC stars in the Tsallis formalism context 
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PACS 97 . 10 . Kc - Stellar rotation 

PACS 97. 10. Yp - Star counts, distribution, and statistics 
PACS 98 . 20 . -d - Stellar clusters and associations 

Abstract - The theoretical distribution function of the projected rotational velocity is derived in 
the context of the Tsallis formalism. The distribution is used to estimate the average (sin i) for a 
stellar sample from the Orion Nebula Cloud (ONC), producing an excellent result when compared 
with observational data. In addition, the value of the parameter q obtained from the distribution 
of observed rotations reinforces the idea that there is a relation between this parameter and the 
age of the cluster. 



Introduction. — The stellar rotation is one of the 
Q most challenging problem in astrophysics because it is a 
j— j complex physical phenomenon. The stellar rotation has 
its origin in the stellar formation, when the angular mo- 
i ^ m entum of the progenitor cloud is transferred to the new 
stars. At this stage, the stars can suffer the influence of 
,— I several physical processes, which can also affect the stellar 



rotation. Such processes include proper orbital motion of 



the cloud, innate differences in internal pressure due to the 
anisotropy of the cloud and the action of the local galac- 
t~--» tic magnetic field on the parent cloud, as well as shock 
waves from supernova explosions. The range of processes 
i-H as discussed makes the initial range of rotation rates diffi- 
04 cult to understand, but when evolution of rotation rates is 
considered along with stellar evolution, processes such as 
magnetic braking and structural changes within the star 
(e.g., development of a radiative core) arc also important. 

5h \ The projected rotational velocity, ysinz, where V is 
, the true (equatorial) rotational velocity and i is the angle 
between the axis of rotation and the line of sight, can be 
derived from the Doppler broadening of spectral lines of 
the star [TH3]. However, it is usually difficult to know 
the inclination of the rotation axis of a single star, being 
possible only in special cases, e.g., when the rotational 
period can be measured (spotted stars). Another special 
case is in the binary systems with short orbital period, in 
which the orbital motion is synchronized with the rotation 
one and rotational and orbital axes are parallel [3J. 

Even given the difficulty in determining the equato- 
rial velocity, one can to derive some global properties of 
the true rotation from a frequency distribution of V sinz. 



However, this requires some assumptions about the fre- 
quency of a given angle of inclination of the rotational 
axis for the group of stars. In general it is assumed that 
the axes of rotation of the stars are randomly oriented (e.g. 
5 8 ). In this context, the relation between the observed 
and true rotational velocity for a sample of stars is always 
(V smi)/(V) = 7r/4 [5], independently of the environment 
of stellar formation. 

The hypothesis of randomly oriented axes was exam- 
ined in several works and, although it has not been ruled 
out, is not always in agreement with the observed velocity 
distributions (see discussion in [TU], § 3.1.2; [TT], § VI). 
In addition, later studies showed a tendency for stars of a 
given spectral type to present rotation correlated with the 
galactic coordinates as well as with the galactocentric dis- 
tance (e.g. [IMI3I)- These studies open up the possibility 
of considering a non-random distribution for the inclina- 
tion of the rotational axes, even if only in special cases. 
In fact, several cases in which a power-law fits better than 
the Gaussian function (e.g. [T3TOfl lTg l lT ^ I3T]), show that 
the main issue may not just be deciding which mathemat- 
ical function can be used but what statistical theory one 
should consider. 

It is now well established that a considerable number 
of stochastic phenomena are better described by power- 
law distributions. For example, one could mention the 
solar neutrino problem [22] . peculiar velocities of galax- 
ies H3] and self-gravitating polytropic systems [23] ■ Fun- 
damentally it could be included in that list of stochastic 
phenomena all the systems whose entropies are more ap- 
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propriately addressed to q-entropy [23], postulated as 




(1) 



where pi denotes the probabilities of the microscopic (in- 
dividual) configurations. In the limit 9=1, the g-entropy 
recovers the standard Boltzmann form. That statistic 
is therefore a generalization of the Boltzmann statistics. 
Usually named noncxtcnsive statistics, the Tsallis thermo- 
statistics states that the entropy of the sum of two systems 
is not simply the sum of individual entropies, but provides 
for cases where such a sum may result in a greater or lesser 
entropy. 

The main point of the considerations presented so far 
is that, besides the exponential distribution considered in 
study as Chandrasekhar & Munch [5] and Deutsch [2j5] , we 
are often faced with stable power-law distributions (the q- 
exponentials in equation [4j , since such distributions are 
ubiquitous in physical phenomena. The ubiquity is due to 
the fact that such distributions naturally obey the general- 
ized central limit theorem (cf. [27]). Thus, Gaussian- type 
diffusion and anomalous diffusion, for example, may be 
unified into a single scenario which is the nonextensivc 
statistics framework [28Tl31j . 

In this paper we present a new way to derive the dis- 
tribution function of the observed stellar rotations in the 
context of the nonextensive statistics. The distribution 
function is tested with a sample of stellar rotation from 
the ONC and accurately reproduces the observational re- 
sults. The paper is organized as follows. In second sec- 
tion we derive the power-law distribution proposed for this 
study and the general equations of its moments. In third 
section we present the work sample and the method to de- 
termine the average (sini) from the observational sample, 
and discuss the results. Fourth section summarizes the 
main results and in the appendix some relevant calcula- 
tions is presented. 

The power-law distribution function. The ro- 
tational kinetic energy of a star with mass M and ra- 
dius R rotating as a rigid body around its own axis is 
E Iot = (l/2)Itu 2 , where I = (2/5)MR 2 is the moment of 
inertia and u = VR^ 1 is the angular velocity of the rotat- 
ing star. Defining (J m - m = (V sini)iJ _1 , it follows that the 
minimal rotational energy, E m i a of the star is given by 
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Consider a group of stars with similar mass, radius and 
age, such that the rotational energy E m - m depends only 
on Vsmi. Suppose further that the energy E m - m is dis- 
tributed according to the function ip(E m i n ). The number 
of stars with energy E m - m between E min and E mnl + d-E m i n 
is 



dN(E n 



p(F 



a } d£min ) 



(3) 



where p{E m \ n ) is the probability density of stars with ro- 
tational energy -E m i n - 



It is well known from observational studies, that the 
probability density of observed rotational velocities van- 
ishes with V sinz. Then it is reasonable to take the distri- 
bution function proposed by Deutsch [32] as an inspiration 
and propose that the probability density decrease accord- 
ing to a generalized exponential law. We are interested 
in deriving the distribution of the rotational velocities in 
the context of the Tsallis formalism. We therefore propose 
that the probability density is governed by the non-linear 
differential equation, dp/dx = ~f3p q , whose solution is 
given by 

p(x) = [l-(l-q)f3x] 1/(1 - q \ (4) 

This function is the (/-deformation of the usual exponential 
function cxp(— fix) and reduces to it in the limit q — > 1 
(cf. |33j). Accordingly, we can rewrite Equation ([3|) as 



dN(E min ) = [1 - (1 - <z)/3£ min ] 1/(1 - 9) dE B 



(5) 



from which it follows that 
2 



dN(y) = -MRy 
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y 



2 1 1/(1-9) 



dy, (6) 



where y = V sin i and <T y = (/3A//5)" 1 / 2 is the half- width 
of the distribution. 

Equation (j5]) gives the number of stars with a given ob- 
served rotational velocity Usini from a sample of stars 
with similar mass, radius and age. It follows, in the for- 
malism proposed by Tsallis, that the distribution function 
of the observed rotation for that group of stars is given, 
apart from a normalization constant, by 

1/(1-9) 



<Pq{v) = V 



l-(l-9)?s 



(7) 



The above equation has also been derived from a different 
way by Soares et al. [H] . 

In accordance with Kraft j34) , we can write the distribu- 
tion function of true rotations, apart from a normalization 
constant, as 

1/(1-9) 



Fg(V) = V 2 



V 2 
v 



(8) 



where ay is the half-width of the distribution. Since the 
function (j^J) recovers the usual exponential in the limit 
q — > 1, Equations Q and jSJ recover the standard form 
of <p(y) and F(V) [261I34] . respectively, when q — > 1. 

Moments of the distributions. The r-th moment of the 
distribution F q (V) is calculated by making 

(V r )= V r F q (V)dV 7 (9) 



where y max = ay/y/1 - q for q < 1, or V max = oo for 
q > 1, in order to ensure the positivity of the function. 
Hence Eq. © provides 



(V r ) = -= 

\ 7T 



ay 
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for q < 1, and 



<n = -7= 



<7y 



7T VV9-1 
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(11) 



for q > 1. The moments of order r = 1, r = 2 and r = 3 
in Equations (fTU|) and (fTTj) give the mean, mean-square 
deviation, and asymmetry of F q (V), respectively. In these 
cases, however, it should be noticed that the convergence 
is achieved only up to certain values of q: q < 3/2 for 
r = 1, q < 7/5 for r = 2 and q < 4/3 for r = 3. 

Similarly to the previous case, one can determine the 
moments of (f q (y) as being 



(y r ) = 



r r 



r 



3-2q 
1-9 



(12) 



and 



(y r > = 



v / 9^1 



-) r (^T 
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for q < 1 and g > 1 respectively. As in the previous case, 
the moments of order r = 1, r = 2 and r = 3 in Equations 
(112[) and (| 1 3|) give the mean, mean-square deviation, and 
asymmetry of <p q (y), respectively. For the regime q > 1, 
the convergence of the moments for y is: q < 5/3 for r = 1, 
q < 3/2 for r = 2, and c/ < 7/5 for r = 3. 

Relation between the moments. Chandrasckhar & 
Munch |J|] had shown that when y = irsini, with (p(y) 
being affected by the randomness factor sini, there is a re- 
lation between the moments of these frequency functions 
which enables us to pass from the moments of the one to 
the moments of other. As is pointed out in the same work, 
the relation is still valid, even when y = xip(i), where 
is an arbitrary integrable function. Using Eqs. (fTU)) , (jTTJ) , 
(|12p and (|13p. we have obtained a similar relation between 
the moments of the distribution of <p q (y) and F q (V): 



(y r ) = ^ 

(V r ) 2 
for q < 1, and 
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for g > 1. The order r = 1, r = 2 and r = 3 in Equa- 
tions dHJ) and ([15]) give the relations between the mean, 
mean-square deviation, and asymmetry of the distribution 
functions F q (V) and <p q (y). It is important to notice that 
in the limit q — > 1 all q-dependent terms in Equation (1151) 



trend to V ( -^j ) and therefore 



(y r ) _V^r(§ + i) 



lim , 



2 r 
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(16) 



where we have assumed a y = cry. Equation ((16]) is 
precisely the relation between moments found by Chan- 
drasekhar & Munch ( 0, Eq. 18). 

The proposed general power-law function (J7J) can de- 
scribe distributions of projected rotational velocities for 
stars from different populations, depending on the value 
of the parameter q. As we will discuss later, the parame- 
ter q constitutes a link between the theoretical distribution 
and the observational one. 

Confronting model with observational data. 

The sample. Our sample was selected out of the data 
provided by Rhode et al. [TD] to match the following cri- 
teria: projected rotation, Vsini, and true rotation, V, 
greater than ~ llkm/s; and sin i < 1. Applying these 
criteria we obtained a sample of 86 stars whose sin i pa- 
rameters are accurately determined. Stellar radii are es- 
timated from combining the luminosities, L, and effective 
temperatures, T e //, displayed in |10) (Table 1), according 
to the relation R oc L 1 / 2 T~ 2 . The stellar parameters L 
and T e ff are uncertain by about 0.2 and 0.02 dex, respec- 
tively (see Hillenbrand [20]). Adding these uncertainties 
in quadrature, we have that the error in R is less than 
0.11 dex. The masses and spectral types for stars in our 
sample can be found in [20] . 

The main features of the stars in our sample are: (i) 
age less than 1-2 Myr; (ii) spectral types ranging from G6 
to M5, with 89% of that stars ranging from K0 to M4; 
(hi) stellar masses covering a range between O.IMq and 
2.7 Mq, with median 0.3 M ; and (iv) stellar radius, R, 
ranging from about 1.2 Rq to 8.2 R Q , with a median of 
2.4Rq. Then our sample consists of stars whose main 
parameters are very similar, as required for this study. 

The projected rotational velocities were obtained from 
the high-dispersion spectra of stars by using the Fourier 
cross-correlation method. The smallest V sin i value that 
is taken to be reliably measured is ll.Okm/s. For a com- 
plete discussion on the observational procedure and error 
analysis of V sin i data, the reader is referred to [TUJ ■ The 
rotations have been calculated from measurements of the 
stellar rotational period, P, and radius, R, using the equa- 
tion V = 2ttRP~ 1 . The calculated rotation values have an 
uncertainty around 11%. The main source of uncertainty 
in these values are the stellar radii, as measured periods 
have a typical accuracy of 1% or better (cf. [TUJ, see also 
Choi & Herbest [T7]). 

The average (sin i) . In this section we show that Func- 
tion ([7]) reproduces accurately the average inclination an- 
gles of the rotational axes, (sini), for stars in the sample. 
First we have determined the values for the parameters 
q and o~ y of the curve that best fits the distribution of 
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the projected rotational velocities. In order to avoid bi- 
ases due to the choice of the bin when using histograms 
to represent frequency distributions, we fit the integral of 
Equation Q to the cumulative distribution of the pro- 
jected rotational velocities. Figure [1] shows the cumu- 
lative distribution and the integral of Equation ([7]) that 
best fits the distribution. We have found q = 1.33 ± 0.03 
and (Ty = 20.9 ± 0.88. To adjust the curve to the dis- 
tribution of V sin i we have used a Levenbcrg-Marquardt 
nonlinear least squares algorithm. The fit parameters are 
X 2 /dof ~ 2.01 x 10~ 3 , for a number of degrees of free- 
dom dof = 84, and probability P > 99%. The region of 
low rotations, by reason of lack of measurements of V sin i, 
presents the largest differences between observational data 
and the model. This lack of the data is due to the im- 
possibility for detecting precise rotational velocities below 
11 km/s in the observational campaign of Rhode et al. [TU] . 
The observed average sine of the inclination of the rota- 
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Fig. 1: Cumulative distribution of the observed rotational data, 
V sini, and the integral curve of the proposed function that 
best fits the distribution. The best-fit parameters are q = 
1.33 ± 0.03 and a = 20.9 ± 0.88. 

tional axes of the ONC stars is (sini) = 0.64 ± 0.07. This 
average was calculated by combining the values of V sin i 
and V for our sample. The uncertainty in this value is 
mainly due to the error of 11%, associated with the rota- 
tion, V. Assuming a random distribution of stellar rota- 
tional axes, the mean value is expected to be n/4 — 0.79, 
which is very different from the observed value. This is 
also displayed in Figure [2j where we represent the distri- 
bution of Vsini as a function of V. In this figure, the 
majority of points are well distributed around the curve 
(Vsini)/(V) = 0.64, but the curve (Vsini)/(V) = tt/4 is 
located near the upper limit of the set of points. Rhode 



et al. [TU] discuss this discrepancy between the observed 
average of sini and theoretical expected value. They con- 
sidered that a possible explanation is that the rotational 
axes in young clusters such as ONC could be aligned rather 
than oriented randomly. In fact, the stars in association 
are formed from a cloud of interstellar material in a rel- 
atively short interval of time. In the formation process, 
the angular momentum of the parent interstellar cloud 
is transferred to the newly-formed stars. It is therefore 
reasonable to assume that the orientation of the angular 
momentum of these stars initially reflect the momentum 
of the rotating progenitor cloud. Thus, a association of 
very young stars tends to have a preferential orientation 
of axes of rotation, although this feature can be lost grad- 
ually over time (cf. [35ll36] h 

i i | i i | i i i | 

q = 1.0 



<sin i> 
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2ttR/P (kms-') 

Fig. 2: Observed rotation, Vsini, as a function of rotation, 
V — 2ttRP~ 1 , for the 86 stars of the sample. The curve q = 1 
represents the function (Vsini)/(V) = 7r/4; (sini) indicates 
(Vsini)/(V) = 0.64; q = 1.33 marks the position of the curve 
(Vsini)/(V) = 0.62; and the last one shows the best-fit curve 
(Vsini)/(V) = 0.51, included here only for comparison with 
the others. The error bars on Vsini are from |10| . 

By using Equation (fT5)) with parameter q = 1.33, as 
obtained from the fit of the V sin i cumulative distri- 
bution, we calculate the expected value for the average 
(sini) = 0.62 ± 0.04. This value is approximately equal to 
the average of sini observed, as shown in Figure [21 where 
the curve (Vsini) /(F) = 0.62 almost coincides with the 
observed one (Vsini)/(V) = 0.64. The expected average, 
calculated in accordance with our model, has a connec- 
tion with the observational data through the parameter q, 
which is derived from the distribution of the projected ro- 
tational velocities. It is therefore understandable that our 
proposal provides an expected average value of (sin i) close 



p-4 



Rotation of ONC stars in Tsallis formalism 



to the observed one. Our model is particularly important 
in cases where only the projected rotational velocities are 
available, because it can provide an average (sini) closer 
the true observed value than can the standard model. The 
standard relation (V suii) / (V) = 7r/4 is appropriate only 
in cases where there is no preferential orientation of stel- 
lar angular momenta. The relation expressed in Equation 
(1151) can be used in any situation, even if the distribution 
is completely random (q = 1). 

In line with the discussion above, if one assumes that the 
younger the cluster the greater the tendency for a preferen- 
tial orientation of stellar angular momenta, the parameter 
q can be associated with the degree of randomness of the 
angular momentum orientations. It can also be related 
to the age of the cluster. Using Equation ([7]), Soares et 
al. |14j have analyzed the distribution of a sample of pro- 
jected rotational velocities of stars in the Pleiades cluster 
(115 Myr) and have found q = 1.38. Utilizing the same 
method, Santoro [21] has found q = 1.51 for the Hyades 
cluster (665 Myr). These results, in combination with the 
present study, reinforce the idea that the parameter q is 
related to the age of the clusters. However, a detailed 
study with stellar rotational data from clusters of different 
ages using the methodology presented here in this work is 
needed in order to determine whether such a relation ex- 
ists. In order to confirm this relation, it is necessary to 
verify that q = 1 for the oldest clusters. 

Conclusions. — We have derived the distribution 
function of the projected rotational velocity in the context 
of the Tsallis formalism. The distribution function is in- 
dependent and recovers the Deutsch's function in the limit 
q — > 1. The parameter q = 1.33 ± 0.03 was determined 
from a sample of 86 stellar projected rotational velocities 
from the ONC. The expect average (sini) = 0.62 ± 0.04 
was calculated using the relation between the first mo- 
ments of the theoretical distribution functions of the ob- 
served, Vsini, and the true, V, rotational velocities in 
the Tsallis formalism. The result reproduces accurately 
the average (sini) = 0.64 ± 0.07 from the observational 
data. The procedure presented in this work constitutes 
an efficient method to determine the relation between the 
moments of the distributions of the observed and the true 
stellar rotational velocities. We suggest and discuss the 
possible existence of a relation between the parameter q 
and the degree of randomness of the angular momentum 
orientations of stars. As well we suggest the possible ex- 
istence of a relation between the parameter q and the age 
of the stellar clusters. However, a detailed study involving 
clusters of different ages is needed in order to determine 
whether any relation exists. This study is being developed 
for a forthcoming paper. 
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